X-ray pulsations from the radio-quiet gamma-ray pulsar in CTA 1 
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ABSTRACT 

Prompted by the Fermi LAT discovery of a radio-quiet gamma-ray pulsar 
inside the CTA 1 supernova remnant, we obtained a 130 ks XMM-Newton ob- 
servation to assess the timing behavior of this pulsar. 

Exploiting both the unprecedented photon harvest and the contemporary 
Fermi LAT timing measurements, a 4.7a single peak pulsation is detected, 
making PSR J0007+7303 the second example, after Geminga, of a radio-quiet 
gamma-ray pulsar also seen to pulsate in X-rays. Phase-resolved spectroscopy 
shows that the off-pulse portion of the light curve is dominated by a power- 
law, non-thermal spectrum, while the X-ray peak emission appears to be mainly 
of thermal origin, probably from a polar cap heated by magnetospheric return 
currents, pointing to a hot spot varying throughout the pulsar rotation. 

Subject headings: stars: neutron — pulsars: individual (PSR J0007+7303) 
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Introduction 



The Fermi LAT discovery (lAbdo et al.l 120081 ) of a pulsed gam ma-ray signal 



the position of the candidate neutron star (NS) RXJ0007.0+7303 (IHalpern et al 



rom 



20041) 



iusid e CTA 1, a 5,000 to 15,000 y old supernova remnant (SNR) at a distance of 1.4±0.3 



kpc ( jPineault et al. 



19931 ). heralded a new era in pulsar astronomy. Besides fostering an 



extraordinary increase in the number of pulsars detected in gamma rays, the Fermi LAT 



first pulsar cat alogue (lAbdo et al. 



as mi 



lisecond ( lAbdo et al. 



2010aj) marks the birth of new pulsar sub -families, su ch 



2009al ) and the radio-quiet gam ma-ray pulsars (lAbdo et al 



2009b|). Indeed, with the recent detection of 8 new objects ( jSaz Parkinson et al. 



2010h the 



gamma-ray dis c overed pulsars now account for one third of the pulsating NSs seen by Fermi 



( jCaraveo 



2010 



Ray fc Saz-Parkinson 



20101 ) . While the detection of so many radio-quiet 
objects points to a gamma-ray beaming covering a solid angle much larger than the radio 
one, one wonders if the difference between radio-quiet and radio-loud gamma-ray pulsars is 
just a geometrical one. 

X-ray astronomy offers an independent way to study NSs by assessing their non-thermal 
magnetospheric emission together with the temperature and emitting area on their surface. 



Previous X-ray studies of CTA 1 centra 



and ROSAT observ ations, 



XMM observations, 



Slane et al 



( l2004h : 



regions unveiled a central filled SNR (ASCA 



Sewardl (119951)) and a poin t sourc e (Chandra and a short 



Halpern et al 



(120041 )). with a jet-like feature, 



embedded in a compact nebula. Standard FFT searches on the XMM data failed to detect 
pulsation, mainly owing to the source faintness. Redoing the exercise using the Fermi 
LAT timing information yielded unconvincing results, thus prompting the request of a long 
XMM-Newton observation. 



1 Based on observations with XMM-Newton, an ESA science mission with instruments 
and contributions directly funded by ESA member states and the USA (NASA). 
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Searching for the source pulsation was, indeed, the main goal of our 130 ks long 
XMM-Newton observation. The unprecedented harvest of X-ray photons, while unveiling 
the Fermi LAT periodicity, allowed also for a new detailed study of the source spectral 
shape. 

Using the newly acquired timing and spectral information, we compare the X-gamma 
behavior of this puzzling radio quiet NS with well known X-ray emitting gamma-ray pulsars. 



2. Observations and data reduction 



Our deep XMM-Newton observation of the CTA 1 system started on 2009, March 7 



at 15:11:10 UT and lasted 130.1 ks. The pn camera ( jStrueder et al. 



20011 ) of the EPIC 



instrument was operated in Small Wind ow mode (time res olution of ~ 5.6 ms over a 4' x 4' 



field of view), while the MOS detectors ( iTurner et al. 



20011 1 were set in Full frame mode (2.6 



s time resolution on a 15' radius field of view). We used the XMM-Newton Science Analysis 
Software v8.0. After standard data processing (using the eppr oc and emproc tasks) and 



De Luca fc Molendi 



200J), 



screening of high particle background time intervals (following 
the good, dead-time corrected exposure time is 66.5 ks for the PN and 93.5 ks for the 
two MOS. The resulting 0.3-10 keV MOS image is shown in Figure [H In order to get 



a sharp view o: 



(IGarmire et al. 



the diffuse emission in the CTA 1 system, we also used a Chandra/ ACIS 



20031 ) observation of the field, performed on 2003, April 13 (50.8 ks observing 
time - such dataset was included in the investigation by Halpern et al. 2004). We retrieved 
"level 2" data from the Chandra Science Archive and used the Chandra Interactive Analysis 
of Observation (CIAO) software v3.2. 
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2.1. Spatial-spectral analysis 



The angular resolution of XMM-Newton telescopes' is not sufficient to resolve the 



pulsar (PSR) from the surrounding pulsar wind ne bula (PWN), T 



Manzali et al 



jus, we used the 



(120071 ) to disentangle 



spatial-spectral deconvolution method developed by 
the point source from the diffuse emission, taking advantage of their different spectra and 
angular distribution. 

1. for each EPIC instrument, we extracted spectra from three concentric regions of 
increasing radii (0-5", 5"- 10", 10"-15"). 

2. based on the well known angular dependence of the EPIC Point Spread Function 
(PSF), we estimated the PSR encircled fraction in each region. Since the target is 
on-axis and most of the counts are below 1 keV, we used PSF model parameters for 
an energy of 0.7 keV and null off axis angle. 

3. We used Chandra data to compute the PWN encircled fraction in each region. To 
this aim, we simulated the PSR point spread function using the ChaRlS and MARx[] 
software packages, assuming the PSR spectrum published by Halpern et al.(2004). 
We positioned the simulated point source at the actual PSR coordinates observed 
by ACIS, and we normalized it in order to match the observed peak counts in the 
ACIS image. Then, we subtracted the simulated image from the observed one. This 
yields an image of the diffuse emission surrounding the pulsar. However, assuming 
pure pulsar emission at the image peak likely results in underestimating the inner 
(r< 1.5") PWN. To correct the residual image for such an effect, we decided to replace 



2 http:/ /cxc. harvard.edu/chart/ 
3 http://space.mit.edu/CXC/MARX/ 
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counts in the inner 1.5" with a poissonian distribution having a mean value equal to 
the average number of counts in the 1.5" — 5" surrounding annulus. As a last step, 
following Manzali et al.(2007), we degraded the angular resolution to match the EPIC 
PSF, obtaining a map of the PWN surface brightness as seen by EPIC. 

4. we fit a two component (PSR+PWN) model to all spectra, freezing the PSR and 
PWN normalization ratios to the results of the previous steps. Uncertainty in best fit 
parameters induced by errors in the encircled fractions is estimated to be negligible 
with respect to statistical errors. 



Accounting at once for both the PSR and the PWN, such an approach yields best 
fit paramet ers for both sp e ctral models. A more detailed description of the method can 



be found in 



Manzali et al. 



( 120071 ). Together with the XMM spectra, we fitted the spectra 
obtained from the pulsar (1.5" circle radius) and the nebula (15" circle radius) in the 
Chandra observation. We used CIAO 4.1.2 software acisspec to generate the spectrum as 
well as the response and effective 

We focus here into the spectral analysis (step 4). Background spectra for each EPIC 
camera were extracted from source-free regions within the same chip. Ad-hoc response 
and effective area files were generated using the SAS tasks rmf gen and arf gen. Since in 
our approach encircled energy fractions for the PSR and PWN are computed a priori and 
then used in the spectral analysis, effective area files are generated with the prescription for 
extended sources, without modelling the PSF distribution of the suorce counts. Spectra 
from the three regions were included in a simultaneous fit using the combination: 



(interstellar absorption) x ( /^(PSR model) + e,(PWN model) ) 
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where pi and are the PSR and PWN encircled fractions within the i extraction 
region. 

The interstellar absorption coefficient does not depend on i. For the PWN, we used a 
power law model. Although the PWN spectrum is expected to vary as a function of the 
position, the relatively small photon statistic prompted us to fit a single photon index T PWN 
to all regions. Of course, PSR parameters do not vary in the different annuli. For the PSR 
emission, we tried a simple power law, the combination of a power law and a blackbody as 
well as the combination of a power law and a magnetized neutron star atmosphere model 
(nsa in XSpec - assuming a 1.4 M neutron star with a radius of 13 km and a surface 
magnetic field of 10 13 G). 

All the models for the PSR emission yield statistically acceptable fits (power law: xt 
=91.5, 124 d.o.f.; blackbody+power law: x^=85.8, 121 d.o.f.; nsa+power law: x^=86.8, 121 
d.o.f.). 

The resulting parameters are summarized in Table 1. As expected, the nsa model 
yields smaller temperatures and larger emitting areas than the blackbody model. 

Using the best fit blackbody+power law model, within a 15" circle in 0.3-10 keV, we 
estimate that 47% of the pn counts come from the PSR, 32% from the PWN and 21% are 
background (instrumental as well as cosmic). 

In order to discriminate between the purely non-thermal and the composite (thermal + 
non-thermal) description of the pulsar emission, the high resolution, temporal information 
provided by the pn instrument is crucial. 
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2.2. Timing Analysis 

4989 pn events in the 0.15-10 keV energy range were extracted fr om a 15" circle, c e ntere d 



Pellizzoni et al 



teoosh . 



on the gamma-ray pulsar. PATTERN selection was performed as by 
X-ray photons' times of arrival were barycentered according to the PSR Chandra position 
(RA 00:07:01.56 , Dec 73:03:08.3) a nd then folded according to an accurate Fermi-LAT 



timing solution ( lAbdo et al. 



2010af ) that overlap our XMM dataset (the pulsar period at 
the start of our XMM-Newton observation is P=0. 3158714977(3) s). Such exercise was 
repeated selecting photons in different energy ranges. 

A 4.7cr pulsation is seen in the 0.15-2 keV energy range (null hypothesis probability of 
1.1 x 10~ 6 , according to a x 2 test), characterized by a single peak, which is out of phase 
with respect to the gamma-ray emission. Light curves computed for different energy ranges 
are shown in Figure El phase aligned with the gamma-ray one. 

Results from the previous section make it possible to compute a net (background and 
PWN - subtracted) pulsed fraction of 85 ± 15% in the 0.15-0.75 keV energy range. No 
pulsation is seen in the 2-10 keV energy range. Assuming a sinusoidal pulse profile, we 
evaluated a 3a upper limit of 57% on the net pulsed fraction. Such a difference in the overall 
source pulsation as a function of the photon energy does not support the single- component 
model for the PSR emission, pointing to the composite thermal plus power law model. 



2.3. Phase-resolved spectroscopy 

Phase-resolved spectral analysis was performed by selecting on- and off-pulse portions 
of the light curve. We selected pn events (15" extraction radius, PATTERN 0) from the 
phase intervals corresponding to the peak and to the minimum of the folded light curve. 
These spectra are plotted in Figure H] where they are seen to differ only at low energy, 
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while they appear superimposed for E> 1.2 keV. We adopted the best fit model computed 
in sect. 12.14 featuring the composite black-body + power-law spectrum for the PS R, as 



a template to descri 



De Luca et al. 



Caraveo et al 



( 120041 ) and 



De the phase-resolved spectra. Following 
( 120051 ). we fixed all spectral parameters (including the PWN component) 
at their phase-averaged best fit values and we used the the normalizations of the PSR 
thermal and non-thermal components to describe the pulse-phase modulation. The spectral 
variation may be well described as a simple modulation of the emitting radius of the 
thermal component, keeping the power law component fixed (x 2 = 26.2, 31 d.o.f. for the 
blackbody+power law model; y 2 = 25.9, 31 d.o.f. for the nsa+power law model). Using the 
blackbody model, the emitting radius of the thermal component varies from 2A2^_\\ 2 m to 



600 



+68 
-75 



m as a function of the star rotation phase. Using the nsa description, the emitting 



radius varies from 560^5gg m to 4380lj^5 m - Quoted values of the emitting radii are as 
measured from a distant observer. Such a variation could easily account for the totality of 
the X-ray pulsation. 

Fitting the on-off spectra using a single power law component does not yield acceptable 
results (x 2 = 48.9, 31 d.o.f.), while the paucity of the counts does not allow to test a model 
where both thermal and non-thermal components are allowed to vary. 



2.4. Extended emission spectral analysis 



Diffuse emis sion, already discovered by ROSAT and ASCA (IHalpern et al. 



Slane et al. 



2004 ; 



20041 ). pervades the entire EPIC/MOS field of view. A thorough analysis of 
such emission, requiring ad-hoc background subtraction/modelling techniques, is beyond 
the scope of this paper, focused the pulsar phenomenology. For completeness, we include 
a simple study of the inner and brighter portion of the diffuse emission (within ~ 150" 
from the pulsar). Using the brightness profile along different radial directions, we selected 
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two elliptical regions (ellipse 1 and 2, see Figure 1) - excluding the inner 15" radius circle 
- and we extracted the corresponding spectra from the MOS data. Background spectra 
were extracted from a region outside ellipse 2. Since ellipse 1 lies within the pn field of 
view, we extracted also a pn spectrum for such a region. Owing to the dimension of the 
pn field of view, the pn background spectrum was extracted from a region within ellipse 2. 
However, the difference in surface brightness between the two ellipses is large enough to 
induce a negligible distortion to the pn ellipse 1 spectrum. The extended emission is well 
described by a power law spectrum with an index of 1.59 ± 0.18 in the inner portion (ellipse 
1, x 2 =123.9, 90 dof) and of 1.80 ± 0.09 in the outer portion (ellipse 2, x 2 =189.3, 137 dof). 
The observed flux is of 1.60±0.09 xl0~ 13 erg cm~ 2 s _1 and of 1.98±0.06 xl0~ 12 erg cm~ 2 
s _1 for the inner and outer portions, respectively. Owing to smaller collecting area as well 
as larger background per unit solid angle, Chandra/ACIS data yield consistent, although 
less constrained, results for such an extended emission. 

3. Discussion and Conclusions 

After Geminga, PSR J0007+7303 now becomes the second example of a radio-quiet 
gamma-ray pulsar also seen to pulsate in X-rays. 

Our deep XMM-Newton observation characterizes the system emission as follows. 

The PWN X-ray spectrum at the position of the pulsar can be described by a power 
law with index T PWN — 1.5 ± 0.3. Diffuse, non-thermal emission with a decreasing surface 
brightness is seen across the EPIC field of view, with a photon index steepening as a 
function of the distance from the pulsar (r = 1.80 ± 0.09 at 1 4- 2.5 arcmin distance). 

The X-ray spectrum of the pulsar is a combination of thermal radiation superimposed 
to a non-thermal power law component. Thermal emission comes from a "hot spot", larger 
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than the polar cap computed for a dipole model for PSR J0007+7303 (about 100 m radius), 
but significantly smaller than the entire surface of any reasonable NS. This is true both 
using the blackbody model and the neutron star atmosphere model (see Table 3), although 
the latter yields a larger emitting area and a lower temperature. Indeed, we warn that 
for some thermally emitting PSR the use of the nsa model turned out to be problematic 
(see e.g. the case of PSR B0656+14 discussed by De Luca et al. 2005). Thus, nsa results 
should be ta ken with cau t ion. A t variance with the majority of X-ray emitting isolated 



pulsars (e.g. 



Kaspi et al. 



( 120061 )). no thermal component from the whole NS surface is 
discernible from the XMM spectrum. Using the blackbody model, the 3a upper limit on 
the temperature of a 10 km radius NS is 5.3 x 10 5 K and 4.8 x 10 5 K using the NSA model. 
This makes PSR J0007+7303 by far the coldest NS for its age interval, suggesting a rapid 
cooling for this young gamma-ray pulsar. 

The detection of X-ray pulsation makes it possible to directly compare the PSR 
J0007+7303 multiwavelength phenomenology with that of other prototypical pulsars. 

With a rotation energy loss of 4.52 xl0 35 er(7s _1 and a kinematic age of 13 ky, 
PSRJ0007+7303 is 50 times younger and 10 times more energetic than Geminga, for 
many years the only known radio-quiet gamma-ray pulsar (for a review, see Bignami 
& Caraveo, 1996). While PSR J0007+7303 is a relatively young, Vela-like pulsar, its 
rotational energy loss is intermediate between Geminga and Vela. Thus, it makes sense 



to compare PSR J0007+7303 wit] 
multiwavelength phenomenology ( 


i Vela and Geminga, two NSs wit' 


i a we 


11 established 


Sanwa 


et al. 


2002: 


Manzali et al. 


2007; 


Caraveo et al. 


2004; 


De Luca et al. 


2005; 


Jackson et al. 


2005 


)• 



Starting from the source flux values, if we consider the ratio between the gamma-ray 
and non-thermal X-ray fluxes, we find a value of (5.6 ±1.1) x 10 3 for CTA 1, to be compared 
with (6.8 ± 0.4) x 10 3 for Geminga and (1.3 ± 0.3) x 10 3 for Vela. Thus, the young and 
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energetic gamma-ray pulsar PSR J0007+7303 is somewhat under- luminous in X-rays, 
joining G eminga and PSR J1 8 36+5925, another radio-quiet pulsar also known as Next 



Geminga (IHalpern et al. 



20071 : 



Abdo et al. 



2010b) 



Turning now to the phase resolved spectral analysis, we note that the peak emission of 
the newly measured single-peak X-ray light curve can be ascribed to a hot spot, apparently 
varying throughout the pulsar rotation. Although the hot spot dimension seems too big to 
be reconciled with the NS polar cap, but too small to account for the entire NS surface, 



the varying thermal contribution is indeed re miniscent of o 



pulsars such as Geminga and PSR B1055-52 (IDe Luca et al 



the b ehavior of middle-aged 



20051 ) . 



Although PSR 0007 is the second radio quiet NS seen to pulsate in X-ray, we note 
that the X-to-gamma-ray observational sequence, which was successfully applied with 
Geminga, has been reversed. Here, for the first time, the basic pulsar properties (P, P) 
were discovered through the gamma-ray emission which made it possible the detection of 
the source X-ray pulsations. 

While X-ray observations of gamma-ray pulsars remain crucial to probe the 
source physics, the superb sensitivity of the Fermi LAT has ma de the discovery of 
periodicities far mor e fruit ful in the gamma-ray band, as shown by 



Abdo et al. 



Saz Parkinson et al 



(boid ) 



(120090 ) and 



This work was supported by contracts ASI-INAF 1/088/06/0 and NASA NIPR 
NNG10PL01I30. 
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Fig. 1. — 0.35-10 keV MOS Imaging. The two MOS exposure-corrected images have been 
added and smoothed with a Gaussian with a Kernel Radius of 13". The two white ellipses 
indicate the extended source and the compact PWN. A 15" blue circle indicates the PSR 
position while two unidentified sources are marked with cyan circles. 
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Fig. 2. — EPIC/pn folded light curves in different energy ranges using photons within a 15" 
radius from the Chandra position. X-ray photons' phases were computed according to an 
accurate Fermi-LAT ephemeris overlapping with the XMM dataset: the pulsar period at the 
start of the XMM observation is P=0. 3158714977(3) s and the P contribution was taken in 
account. PATTERN events have been selected in the 0.15-0.35 keV energy range while 
PATTERN <4 have been used in the 0.35-2 keV range. The upper panel shows the LAT 



light curve of the CTA 1 pulsar from 
have been aligned in phase. 



Abdo et al. 



fcoiOal ) to which the XMM light curves 
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data and folded model 




Energy £keV) Energy (keV) 

Fig. 3.— PN and MOS spectra of PSR J0007+7303. PATTERN PN events and PATTERN 
<12 MOS events have been selected among photons within 15" from the target position. The 
spectra are rebinned in order to have at least 25 counts per bin and no more than 3 spectral 
bins per energy resolution interval. The black and the red curves show respectively the PN 
and MOS data and spectral fits. Cyan square-marked curve shows the blackbody component, 
while pulsar power law is shown with blue curves and PWN one with green asterisks. 
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Phase 

0.5 1.0 1.5 




Energy (keV) 

Fig. 4.— Upper panel: X-ray folded light curve of PSR J0007+7303 (0.15-2 keV). Lower 
panel: X-ray spectra relative to the phase intervals shaded in gray (on-pulse, black line) and 
red (off-pulse, red line). Both spectra were fitted with a 3-component model to account for 
the NS thermal emission (cyan symbols) and power law (blue dotted line) as well as the 
PWN power law (green dotted line). With the power law contributions unchanged in the 
two spectra, a significant thermal component is present only in the on-pulse spectrum (cyan 
triangles) while it appears suppressed in the off-pulse one (cyan squares). 



- 19 - 



Table 1: CTA 1 Pulsar and nebula spectra. 





Pulsar 


J0007+7303 




Inner PWN 


Outer PWN 


Parameter 


PL 


PL+BB 


PL+NSA 






N„(10 21 ) 
Tpwn 

FpSR 

kT(keV) 
r L4fcpc (km) 

x 2 


0.63^;1 

i ocr+0.17 
i - zo -0.15 

1 ^fi+ - 16 
i -' 3O -0.14 

- 

91.56 


i-66±8:S 

i rq+0.33 

1.30±0.18 

n i Qo+0.032 

64+ - 88 

u - u ^-0.20 

85.81 


i co+1.19 
i -°°-0.73 

1 4q+ - 32 

i -^ y -0.24 
i 9 c-+0.20 

0-054±8;g» 
86.82 


i.67ii;| 

1.59±0.18 
- 
- 

123.89 


1 oq+0.30 
l.OO_ 27 

1.80±0.09 
- 
- 

189.30 


d.o.f. 


124 


121 


121 


90 


137 


Total Flux . 3 -iofceV a 


12.00±0.10 


13.90±0.36 


14.01±0.41 


16.0±0.9 


198±6 


Total Flux 2 _iofcev a 


O OO+0.37 

S.°o_ 028 


8 69 +0 - 97 

°- uy -0.86 


8.74l ;° 2 


10.1±0.6 


105±5 


PSR Fluxo.3_iofceV a 


6.54±0.53 


8.41±0.98 


8.41±1.00 






PSR FluX2_l0fceV ° 


3.98±0.72 


4 30+ 1 - 62 

^• ou -0.61 


4 32 +1 - 67 

^• OZ -0.65 






Thermal Flux a 




1.55±1.01 


1.68±1.11 







Note. — X-ray spectrum of the pulsar and the nebula. Inner and Outer PWN correspond to emission 
from ellipse 1 and ellipse 2, respectively (see taxt and Figure 1). For the pulsar we provide the power law, 
the blackbody + power law and the magnetized neutron star atmosphere model (nsa) + power law spectral 
fits. Temperatures and emitting radii are as measured from a distant observer. For the pulsar and inner 
nebula we used Chandra, MOS and PN data, while for the outer PWN spectra we used only data from 
MOS1+2 instruments owing to the small FOV of PN. 
a Fluxes are in units of 10 -14 erg/cm 2 s 



